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Probing the Reheating Temperature at Colliders and with Primordial Nucleosynthesis 
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Considering gravitino dark matter scenarios with a long-lived charged slepton, we show that 
collider measurements of the slepton mass and its lifetime can probe not only the gravitino mass 
but also the post-inflationary reheating temperature Tr. In a model independent way, we derive 
upper limits on Tr and discuss them in light of the constraints from the primordial catalysis of 
6 Li through bound-state effects. In the collider- friendly region of slepton masses below 1 TeV, the 
obtained conservative estimate of the maximum reheating temperature is about Tr = 3 x 10 9 GeV for 
the limiting case of a small gluino-slepton mass splitting and about Tr = 10 s GeV for the case that 
is typical for universal soft supersymmetry breaking parameters at the scale of grand unification. 
We find that a determination of the gluino-slepton mass ratio at the Large Hadron Collider will test 
the possibility of Tr > 10 9 GeV and thereby the viability of thermal leptogenesis with hierarchical 
heavy right-handed Majorana neutrinos. 

PACS numbers: 98.80.Cq, 95.35.+d, 12.60.Jv, 95.30.Cq 
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INTRODUCTION 

Big bang nucleosynthesis (BBN) is a powerful tool to 
test physics beyond the Standard Model; cf. [1-7] and 
references therein. Indeed, in supersymmetric (SUSY) 
theories, severe constraints appear owing to the existence 
of the gravitino G which is the gauge field of local SUSY 
transformations and whose mass is governed by the SUSY 
breaking scale. As the spin-3/2 superpartner of the gravi- 
ton, the gravitino is an extremely weakly interacting par- 
ticle with couplings suppressed by inverse powers of the 
(reduced) Planck scale M P = 2.4 x 10 18 GeV [8]. Accord- 
ingly, once produced in thermal scattering of particles in 
the hot primordial plasma [9-15, 17], unstable gravitinos 



5 TeV have long lifetimes, > 100 s, 



with a mass rr. 
and decay during or after BBN. Since the decay prod- 
ucts affect the abundances of the primordial light ele- 
ments, successful BBN predictions imply a bound on the 
reheating temperature after inflation Tr which governs 
the abundance of gravitinos before their decay [9-17]: 
T R < 10 s GeV for m d < 5 TeV [5, 7]. 

We consider SUSY extensions of the Standard Model 
in which the gravitino G is the lightest supersymmet- 
ric particle (LSP) and a charged slepton l\ — such as 
the lighter stau t\ — the next-to-lightest supersymmet- 
ric particle (NLSP). Assuming R-parity conservation, the 
gravitino LSP is stable and a promising candidate for 
dark matter; cf. [10, 12-23] and references therein. Be- 
cause of the extremely weak interactions of the grav- 
itino, the NLSP typically has a long lifetime before it 
decays into the gravitino. If these decays occur dur- 
ing or after BBN, the Standard Model particles emitted 
in addition to the gravitino can affect the abundances 
of the primordial light elements. For the charged slep- 
ton NLSP case, the BBN constraints associated with 
hadronic/electromagnetic energy injection have been es- 
timated [19-22]. Taking into account additional con- 
straints from large-scale-structure formation, apparently 



viable gravitino dark matter scenarios had been identi- 
fied (see, e.g., benchmark scenarios A\p-Ci i m Ref- [22]) 
which are attractive for two reasons: (i) Tr > 10 9 GeV 
is possible so that thermal leptogenesis with hierarchical 
right-handed neutrinos [24-28] 1 remains a viable expla- 
nation of the baryon asymmetry [12, 14, 21, 22, 36] and 
(ii) the gravitino mass can be close to the NLSP mass 



i.e., 0.1 mj < 



so that a kinematical 
determination appears viable [37-39]. With a kine- 
matically determined tHq, one would be able to measure 
the Planck scale Mp at colliders [37-39] and to test the 
viability of thermal leptogenesis in the laboratory [14]. 
Indeed, an agreement of the Mp value determined in col- 
lider experiments with the one inferred from Newton's 
constant [40] G N = 6.709 x 10" 39 GeV~ 2 would provide 
evidence for the existence of supergravity (SUGRA) in 
nature [37]. 

Already in the early paper [41] it had been realized that 
long-lived negatively charged massive particles X~ (such 
as long-lived Z^'s) can form primordial bound states, 
which can affect the abundances of the primordial light 
elements. It was however only recently [42] when it was 
realized that bound-state formation of X~ with 4 He can 
lead to a substantial overproduction of primordial 6 Li via 
the catalyzed BBN (CBBN) reaction 



( 4 HeY- 



D 



3 Li + X- 



(1) 



With an X abundance that is typical for an electrically 
charged massive thermal relic [18], this reaction becomes 



1 Note that flavor effects [27-30] do not change the lower bound 
Tr > 10 9 GeV required by successful thermal leptogenesis with 
hierarchical right-handed neutrinos [27, 28]. However, in the case 
of (nearly) mass-degenerate heavy right-handed Majorana neu- 
trinos, resonant leptogenesis can explain the baryon asymmetry 
at smaller values of Tr [31-34]. Another example for a frame- 
work in which the limit Tr > 10 9 GeV is relaxed is non-thermal 
leptogenesis; see, e.g., [35] and references therein. 
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so efficient that an X~ lifetime of t x - > 5 x 10 3 s is 
excluded by observationally inferred values of the pri- 
mordial 6 Li abundance [42]. In the considered gravitino 
LSP scenarios, this bound applies directly to the lifetime 
of the h NLSP [7, 15, 23, 42-48]: % < 5 x 10 3 s. 2 This 
implies Trig < 0.1 mj in the collider-friendly mass range 
of mj < 1 TeV. Accordingly, the region 0.1 nij < niQ < 
mj , in which the kinematical mg determination appears 
feasible, seems to be excluded by BBN constraints [44]. 
Moreover, within the framework of the constrained mini- 
mal supersymmetric Standard Model (CMSSM), we have 
found that the 6 Li constraint implies the upper limit 
Tr < 10 7 GeV in gravitino dark matter scenarios with 
unbroken R-parity and typical thermal Ti NLSP relic 
abundances [15, 47, 48]. For a standard cosmological 
history, this finding clearly disfavors successful thermal 
leptogenesis within the CMSSM in the case of hierarchi- 
cal heavy right-handed Major ana neutrinos. For a grav- 
itino LSP mass range that is natural for gravity-mediated 
SUSY breaking, the 6 Li constraint can even point to a 
CMSSM mass spectrum which will be difficult to probe 
at the Large Hadron Collider (LHC) [15, 43, 47, 48]. 

This letter provides a model independent study of 
gravitino LSP scenarios with a charged slepton NLSP 
that has a 6 Li-friendly lifetime of < 10 4 s and a 
collider-friendly mass of mj < 1 TeV. In contrast 
to [15, 47, 48, 55], the investigation presented in this 
work is not restricted to a constrained framework such 
as the CMSSM or to the gravitino mass range suggested 
by gravity-mediated SUSY breaking. Thereby, generic 
results are obtained with a range of validity that includes 
models with gauge-mediated SUSY breaking and/or non- 
standard mass spectra. 

The remainder of this letter is organized as follows. In 
the next section we review that rag could be determined 
by measuring mj and tj- at colliders. We give also an up- 
per limit on mg that depends on and . This limit 
allows us to derive a lower limit for the gravitino density 
fig. Since f2g cannot exceed the dark matter density 
Odm, this leads to conservative Tr limits which can be 
probed in measurements of mj and tj (and the gluino 
mass nig). We use the derived expressions to translate 
the tj- constraint from CBBN of 6 Li into robust upper 
limits on Tr that depend on the mass ratio nig / mj . 
Finally, we show that the requirement Tr > 10 9 GeV 
needed for successful standard thermal leptogenesis pro- 
vides an upper limit on this ratio and thereby a testable 
prediction for LHC phenomenology. 



2 While numerous other bound-state effects can affect the abun- 
dances of 6 Li and the other primordial light elements [43, 49-53], 
the approximate Tr bound is found to be quite robust; cf. [54]. 
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FIG. 1: Contours of m~ n (solid lines) in the plane spanned 
by tt i and mr . BBN constraints from D and 6 Li are indicated 
by the short-dash-dotted (blue in the web version) and by the 
long-dash-dotted (red in the web version) lines, respectively. 
The 6 Li constraints are obtained from the CBBN treatment 
of [47] for upper limits on the primordial 6 Li/H abundance of 
2 x KT 11 (left line) and 6 x HT 11 (right line). 

PROBING m 5 AT COLLIDERS 

In the considered SUSY scenarios with the gravitino 
LSP being stable due to R-parity conservation, 3 the 
charged slepton NLSP has a lifetime tj- that is governed 

by the decay l\ —> Gt and thus given by the SUGRA 
prediction for the associated partial width 

487rm 2 ~M 2 / mf=\ "* 487rm?-M 2 

n k— - 1 - Hm > 2-^ > ( 2 ) 

h \ hJ h 

where the rightmost term underestimates tj by at most 
5% (30%) for mg < O.lmj; (m s < 0.25 m~). Accord- 
ingly, using Mp = 2.4 x 10 18 GeV as inferred from New- 
ton's constant, the gravitino mass mg can be determined 
by measuring both mj and tj at future colliders [57]. 
Moreover, from the rightmost expression in (2), one can 
extract an upper limit for the gravitino mass 

m d < 0.41 GeV (^)\^)^ ^ , (3) 
which turns into an equality for mg <C mj . 



3 For the case of broken R-parity, see, e.g., [56]. 
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In Fig. 1 contours of m-" 3 ^ (solid lines) are shown in 
the plane spanned by Tr and mr . Only mr > 100 GeV is 

*> 1 » 1 * 1 

considered since long-lived charged sleptons should have 
otherwise been observed already at the Large Electron- 
Positron Collider (LEP) [40]. 

The potential for collider measurements of mj^ is 
promising since each heavier superpartner produced will 
cascade down to the l\ NLSP which will appear as a 
(quasi-) stable muon-like particle in the detector [58-60] . 
For slow Zi's, the associated highly ionizing tracks and 
time-of-fiight measurements can allow for a distinction 
from muons [57-60]. With measurements of the l\ ve- 
locity (3j = vj /c and its momentum pj = \pj |, mj 

can be determined: = p^(l - 0~ ) 1 ^ 2 //% 1 [57]. For 

the upcoming LHC experiments, studies of hypothetical 

scenarios with long-lived charged particles are actively 

pursued [61-63]. In Ref. [61], for example, it is shown 

that one should be able to measure the mass mr of a 

1 

(quasi-) stable slepton quite accurately at the LHC. 

The experimental determination of tj- will be substan- 
tially more difficult than the measurement. If some 
of the sleptons decay already in the collider detectors, 
the statistical method proposed in [57] could allow one to 
measure tj- . Moreover, ways to stop and collect charged 
long-lived particles for an analysis of their decays have 
been proposed for the LHC and for the International 
Linear Collider (ILC) [37-39, 64-69]. These challenging 
proposals could lead to a precise measurement of tj- . 

In addition, these proposals could help to distinguish 
the case of the gravitino LSP from the one of the axino 
LSP [39, 70], which is — as the fermionic superpartner 
of the axion another well-motivated dark matter candi- 
date [23, 71] . Indeed, also the axino LSP can be produced 
thermally in the early Universe [72, 73] and can be as- 
sociated with a h NLSP with O(ms) < < O(days) 
[70, 74]. In the axino LSP case, however, measurements 
of tj- and mr will probe the Peccei-Quinn scale f a [70] 
instead of mg. Keeping in mind that the axino LSP could 
mimic the gravitino LSP at colliders (at first sight), we 
assume in the remainder of this work that it is the grav- 
itino LSP scenario that is realized in nature. 



to their extremely weak interactions, 4 the thermally pro- 
duced gravitino density rjT P depends basically linearly on 
Tr — cf. (5) below — and thus can serve as a thermome- 
ter of the earliest moments of the radiation-dominated 
epoch [10, 12-15, 17]. In particular, this allows for the 
derivation of upper limits on T R since the relic gravitino 
density is bounded from above by the observed dark mat- 
ter density [40]: 



nfh 2 < n d h 2 < n DM h 2 



0.1, 



(4) 



where h ~ 0.7 denotes the Hubble constant in units of 
100 kmMpc _1 s _1 . 

Aiming at a lower limit of fig to arrive at a truly con- 
servative upper limit on Tr, we do not take into account 
the model dependent contributions from decays of scalar 
fields such as the inflaton [75-77] or the ones from NLSP 
decays [18, 20, 78], 5 which can be negligible anyhow, in 
particular, for rag < 1 GeV. Indeed, we focus on the 
SUSY QCD contribution to the thermally produced grav- 
itino density, as derived in the gauge- invariant calculation 
of [13, 14], 



"<3 h SU(3) C 



= 0.117 5s 2 (T R ) 



1 + 



m?(T R ) 



3m| 



In 



1.271 



9s(Tn) 



( T R 



100 GeV/ V!0 10 GeV 



(5) 



with the strong coupling g s and the gluino mass 
nig to be evaluated at the scale given by Tr, 
i.e., g s (T R ) = [g- 2 (M z ) + 3HT r /M z )/(8tt 2 )]- 1 / 2 
and mg(T R ) = [g s (T R )/g s (M z )} 2 mg{M z ), where 
g 2 {M z )/{ATi) = 0.1172 at M z = 91.188 GeV. 

By discarding the electroweak contributions, f2~ p is 
underestimated typically by 20%-50% depending on the 
size of the gaugino masses in the electroweak sector [14- 
16]. Thus, these additional contributions would tighten 
the Tr limits derived below only by a factor that is not 
smaller than about 2/3. 

For high-temperatures, 10 6 GeV < T R < 10 10 GeV, 
where (5) derived in the weak coupling limit g s <C 1 is 
most reliable [13, 14], we do now derive a lower limit 



nf»h 2 <nfh 2 \ slJ{3)c 

by manipulating expression (5) as follows: 



(6) 



(i) We use the replacement 



PROBING T R AT COLLIDERS 



In 



1.271 



9s (Tr) 



In 



1.271 



5s (10 6 GeV) 



= 0.256 . (7) 



Within a standard cosmological history, the gravitino 
LSP can be produced in decays of scalar fields such as 
the inflaton [75-77] , in thermal scattering of particles in 
the primordial plasma [10, 12-15, 17], and in NLSP de- 
cays [18-20, 78]. Focussing on the case mg ;§> 100 eV 
in which gravitinos are never in thermal equilibrium due 



4 In gauge-mediated SUSY breaking scenarios, light gravitinos can 
be viable thermal relics if their abundance is diluted by entropy 
production, which can result, for example, from decays of mes- 
senger fields [79—83]. 

J For Tr limits obtained by taking into account contributions to 
U d from NLSP decays, see e.g. [14, 15, 18, 21, 22, 36, 48, 84]. 
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Thereby, 0~ p /i 2 |su(3) IS underestimated at most 

Cr *> 'C 

at Tr = 10 10 GeV and there by a factor of about 
5/8. 

(ii) We use 5 S (10 10 GeV) = 0.85 in the numerator of (5) 
and to evolve m,g(T R ) to the weak scale. Thereby, 
f2~ p /i 2 |su( 3 ) is underestimated at most at Tr = 
10 6 GeV and there by a factor of about 2/5. 

(hi) We use the constant c > 1 to parametrize vtig(Mz) 
in terms of the mass of the slepton NLSP: 



m-g{M z ) — cmr 



(8) 



For concreteness, we use rrig(Mz) to represent the 
gluino mass at the weak scale. While the running 
mass m-g decreases when evolved to higher energy 
scales, we assume implicitly m-g > (even in 
the limiting case c ~ 1) at least up to energy 
scales accessible at the LHC and thereby up to 
temperatures^ well above the freeze-out tempera- 
ture of the h NLSP. For example, since nig de- 
creases by about 20% when evolved from Mz to 
10 TeV, ^~ P /i 2 |su(3) c can thereby be underesti- 
mated for c ~ 1 by a factor of about (4/5) 2 . 

(iv) We drop the contributions from the spin 3/2 com- 
ponents of the gravitino which are given by the 
term in (5) that is independent of m g . Focussing 
on mg < 0.1 m^, the relative importance of the 
spin 1/2 components is minimal for mg ~ 0.1 m-^, 
c ~ 1, and Tr = 10 10 GeV, where the associated 
second term in the first bracket of (5) becomes 
about 8 so that S!g P /i 2 |su(3) c is underestimated 
without the spin 3/2 term by a factor that is not 
smaller than about 8/9. 

Accordingly, we find a guaranteed gravitino density of 



0.174 



1 GeV 



100 GeV J V 10 10 GeV 



T, 



R 



For mg < 0.1 m Ti , 10 8 GeV (10 9 GeV) < T R < 10 10 GeV 
and cmj > 200 GeV (100 GeV), this expression is 

associated with 1.57 < ^| P |su (3 ) c /^g in < 2.1 (1.9). 
Accordingly, the limits derived below may be consid- 
ered to be 'too relaxed' (or 'too conservative') by at 
least a factor of about 1.5. Indeed, 0~ P | SU (3) /^g m 
increases by decreasing Tr and/or cra^ or by increas- 
ing mg so that the corresponding limits presented below 
will become even more conservative. For example, for 
T R = 10 6 GeV, cm Ti ~ lOOGeV, and mg < 0.05 m ri , 
0~ p | SU (3) c /0~ m — 2.4 is encountered so that the asso- 
ciated limits derived below may be considered as 'too 
relaxed' by at least a factor of about 2.4. 



With m~ ax from (3), expression (9) leads to a lower 
limit on Op; in terms of mr and Tr : 

Ct tl tl 



Og/i 2 > 0.422 c 2 



10 4 s\Y 100 GeV V/ T R 



mr 



,10 10 Gcvy ' 

(10) 

Comparing this lower limit with the dark matter con- 
straint (4), we arrive at a conservative Tr limit that can 
be determined at colliders by measuring the slepton mass 
, its lifetime tj i , and c (or the gluino mass m-g) : 



Tr < 



2.37 x 10 9 GeV /On M /i 2 



0.1 



This is one of the main results of this letter. Note that 
this limit can be refined easily. Once an additional con- 
tribution O x to Odm — such as an axion density or a non- 
thermally produced gravitino density ~ TP — is taken for 
granted, the resulting tighter Tr limits can be obtained 
from (11) after the replacement: Odm — > Odm — O x . In 



fact, based on the interplay between 0~ p and the contri- 
bution from NLSP decays, 0~ TP -dependent upper limits 

Cr 

on Tr have been derived to be tested at colliders [84]. 

At this point, one has to clarify to which Tr definition 
the limit (11) applies. The analytic expression (5) is de- 
rived by assuming a radiation-dominated epoch with an 
initial temperature of Tr [13, 14]. In a numerical treat- 
ment, the epoch in which the coherent oscillations of the 
inflaton field dominate the energy density of the Uni- 
verse can also be taken into account, where one usually 
defines Tr in terms of the decay width of the infla- 
ton field [4, 15]. In fact, the numerical result for 0~ p 
has been found to agree with the corresponding analytic 
expression for [15] 



90 


■'V 


U(7rK 2 _ 





1.8 



(12) 



(9) which satisfies ~ 1.8H Ta d(T R ) with the Hubble pa- 



rameter iT r ad(T) = y/g*(T)iv 2 /90T 2 /M P and an effec- 
tive number of relativistic degrees of freedom of g* (Tr) = 
228.75. Thus, (12) provides the Tr definition to which 
the upper limit (11) applies. For an alternative Tr defi- 
nition given by T = £H rad (T^ ] ), 



90 



1/4 



/Mr 



(13) 



the upper limit (11) can be translated accordingly 
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FIG. 2: Contours of Tr ax imposed by Q d h 2 < Sl DM /i 2 < 
0.126 for c = 1 (solid lines) and c = 7 (dashed lines) in the 
plane spanned by tj- and mj . The BBN constraints are iden- 
tical to the ones shown in Fig. 1. 

In particular, the associated numerically obtained f2~ p 
can be reproduced with the analytical expression after 
substituting T R with y/^/l.8T^ ] [15]. 

In Fig. 2 contours of the upper limit Tp nax given in (11) 
are shown for c = 1 (solid lines) and c — 7 (dashed lines) 
as obtained with f^DM^ 2 < 0.126. 6 The most conserva- 
tive upper limit is represented by the limiting case 7 c = 1 
which holds even without insights into rnj. Indeed, in 
this exceptional case, we assume implicitly a mass differ- 
ence, Am = trig — to^ > 0, such that l\-g coannihilation 
effects are negligible and such that long-lived gluinos do 
not appear (avoiding thereby possibly severe hadronic 
BBN constraints from late decaying <?'s [36]). 

With experimental insights into mj i and rrig (and 
thereby into c), the limit (11) can become considerably 
more severe than in the c = 1 case. In Fig. 2 this is 
illustrated for c = 7 which is a value that appears typ- 
ically in constrained scenarios with universal soft SUSY 
breaking parameters at the scale of grand unification 
M GUT ~2x 10 16 GeV. For example, in the CMSSM 
in which the gaugino masses, the scalar masses, and the 



6 This Q,-£,y[h? value is used to allow for a convenient compari- 
son with our previous works [14, 15]. There Qj^j^h 2 < 0.126 is 
understood as a nominal 3<r limit derived with a restrictive six- 
parameter "vanilla" model from the three year data set of the 
Wilkinson Microwave Anisotropy Probe (WMAP) satellite [85]. 

7 This limiting case has already been considered in the benchmark 
scenarios Ci,2 m ^ e f- [22] which are now understood to be dis- 
favored by the 6 Li constraint from CBBN [42]. 



trilinear scalar couplings are assumed to take on the re- 
spective universal values mi«, mo, and Aq at Mqut, 
the ti NLSP region has been found to be associated with 
m~ < 0.21 j2 and thus with c > 6 over the entire 
natural parameter range [47]. 

The Tr limits derived above are in spirit similar to the 
mg- dependent Tr, limits given in Refs. [10, 15, 17, 18, 
21, 22, 84]. Also the sensitivity of these Tr limits on 
the gaugino masses has already been discussed and used 
to provide mg-dependent upper limits on the gaugino 
masses for given values of Tr [12-14, 36]. In fact, once 
the gaugino masses and the gravitino mass m^ (inferred 
from mj and rr) are known, numerical results such as 
the ones shown in Fig. 2 of Ref. [15] will provide a Tr limit 
that is more restrictive than our analytic estimate (11). 
However, as long as insights into a SUSY model possibly 
realized in nature are missing, we find it important to 
provide a conservative and robust estimate of T R nax that 
is insensitive to details in the SUSY spectrum (other than 
the assumption of the G LSP and a l\ NLSP). In partic- 
ular, the comparison of (11) with more model dependent 
limits— obtained, e.g., within the CMSSM [15, 47, 48]— 
demonstrates very clearly the impact of restrictive as- 
sumptions on the soft SUSY breaking sector. Moreover, 
the derivation of T^ ax contours as a function of and 
Tf may turn out to become very useful since mj and 
are the quantities that might be directly accessible 
in collider experiments. Finally, as demonstrated below, 
the presentation of 7 1 ™ ax in the plane spanned by tt and 
mj allows also for a convenient analysis of 7 1 ™ ax in light 
of recent BBN constraints. 

Before proceeding, we would like to stress that the re- 
lation of T™ ax to quantities mj , Tr , and c (or m-g), which 
could be accessible at future collider experiments, relies 
crucially on assumptions on the cosmological history and 
the evolution of physical parameters. For example, for 
a non-standard thermal history with late-time entropy 
production, the thermally produced gravitino abundance 
can be diluted f2~ p — > fi~ p /(S by a factor S > 1 so 
that Tg** -> 5T™ ax [15]. Moreover, the limit (11) can 
be evaded if the strong coupling g s levels off in a non- 
standard way at high temperatures [86] . This emphasizes 
that (11) relies on the assumptions of a standard cosmo- 
logical history and a strong gauge coupling that behaves 
at high temperatures as described by the renormaliza- 
tion group equation in the minimal supersymmetric Stan- 
dard Model (MSSM), i.e., as described by g s (T R ) given 
below (5) . While tests of these assumptions seem inac- 
cessible to terrestrial accelerator experiments, futuristic 
space-based gravitational-wave detectors such as the Big 
Bang Observer (BBO) or the Deci-hertz Interferometer 
Gravitational Wave Observatory (DECIGO) [87] could 
allow for tests of the thermal history after inflation and 
could even probe Tr [88, 89] in a way that is complemen- 
tary to the approach presented in this letter. 
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PROBING T R WITH PRIMORDIAL 
NUCLEOSYNTHESIS 

Gravitino LSP scenarios with a long-lived charged slep- 
ton NLSP can affect BBN as described in the Introduc- 
tion. Particularly severe is the catalytic effect of ( 4 HeZj~)- 
bound states on the primordial abundance of 6 Li [42]. 
Indeed, the CBBN reaction (1) can become very efficient 
at temperatures T ~ 10 keV depending on the l\ abun- 
dance at that time. Observationally inferred upper lim- 
its on the primordial 6 Li/H abundance 6 Li/H| p (cf. [90]) 
can thus be translated into tj- -dependent upper limits on 

the thermal relic abundance of the negatively charged Zj" 
NLSP [7, 42, 45-47, 91]. For a given SUSY model, this 
abundance can be calculated, for example, with the com- 
puter program micrOMEGAs 2.1 [92]. By confronting the 
obtained abundance with the tj- -dependent upper limits 
(see, e.g., Fig. 1 in Ref. [47]), one can extract an upper 
limit on tj [7, 15, 42-47, 91]. This tj limit can then be 
used to determine T™ ax directly from (11). 

Let us perform this procedure explicitly for the follow- 
ing thermal relic l\ NLSP abundance after decoupling 
and prior to decay [18, 36]: 

nr ., / mr \ 

Yr =J±=2Y T -= 0.7 x 10~ 13 — £ - , (15) 

where s denotes the entropy density and the total 

1 1 number density assuming an equal number density of 
positively and negatively charged l\'s. 

We confront (15) with the Yj- limits that emerge 
from a calculation of the 6 Li abundance from CBBN, 
6 Li/H|cBBN, which uses the state-of-the-art result of the 
catalyzed 6 Li production cross section [45] and the Boltz- 
mann equation (instead of the Saha type approximation) 
to describe the time evolution of the ( 4 HeZj~)-bound-state 
abundance [47, 91]. 8 Contour lines of 6 L1/H|cbbn ob- 
tained in this calculation are shown, e.g., in Fig. 1 of 
Rcf. [47]. 9 Working with an upper limit on the primor- 
dial 6 Li/H abundance of [93] 

6 Li/H| p < 2 x HT 11 , (16) 

the corresponding contour given in that figure is the one 
that provides the upper limit for Yr- . As a second more 
conservative limit, we consider in this work also 

6 Li/H| p < 6 x HT 11 , (17) 



8 We thank Josef Pradler for providing us with the 6 Li/H|cBBN 
data from the CBBN treatment of Ref. [47]. 

9 Comparisons of [47] with [7, 43, 46, 53], in which also the possible 
destruction of 6 Li due to li decays is considered, show that those 
effects affect the 6 Li constraint only marginally. 



which provides a more relaxed Yr- limit. 

'i 

In Figs. 1 and 2, the long-dash-dotted (red in the 
web version) lines show the constraints obtained by con- 
fronting (15) with the Yj- limits for 6 Li/H| p < 2 x 10~ n 

(left line) and 6 x 10~ n (right line). It is the region 
to the right of the respective line that is disfavored 
by 6 Li/H|cBBN exceeding the limit inferred from ob- 
servations. Only the constraint from the primordial D 
abundance on hadronic energy release [4, 6] in l\ de- 
cays [20-22] can be more severe than the CBBN con- 
straints [7, 15, 23, 43, 44, 46]. 10 In Figs. 1 and 2, the 
associated disfavored region is the one above the short- 
dash-dotted (blue in the web version) line. For details on 
this constraint, see [15, 22] and references therein. 

In Fig. 1 one sees explicitly that the BBN constraints 
disfavor the region 0.1 < mg < as already em- 
phasized in the Introduction. Moreover, the 6 Li con- 
straint imposes for mg > 10 GeV the lower limit mj > 
400 GeV which implies rrig(M z ) > 2.4 TeV for c > 6. 
This is in agreement with what had previously been re- 
alized in the CMSSM (where c > 6 in the n NLSP re- 
gion [47]): For mg in the range that is natural for gravity- 
mediated SUSY breaking, the 6 Li constraint can point to 
a SUSY mass spectrum that will be difficult to probe at 
the LHC [15, 43, 47, 48]. On the other hand, Fig. 1 shows 
very clearly that the 6 Li constraint becomes negligible for 
rng < 1 GeV, i.e., for the mg range that is natural for 
gauge-mediated SUSY breaking scenarios. 

Figure 2 allows us to read off the T^ ax value imposed 
by the 6 Li constraint in the collider-friendly region of 
< 1 TeV. Our most conservative and thereby most 
robust limit is the one obtained in the limiting case c = 1 . 
There one finds the T£ lax = 3 x 10 9 GeV contour in the 
region that is allowed by the BBN constraints. If the 
gluino turns out to be significantly heavier than the l\ 
NLSP, the T^ ax value will become considerably more se- 
vere. For example, for the case of c = 7, shown by the 
dashed lines in Fig. 2, one finds that a reheating tem- 
perature above T^ ax = 10 8 GeV is disfavored by BBN 
constraints. Note that these limit are quite robust since 
they emerge from the conservative limit (11). In fact, 
one may consider these limits as overly conservative by 
at least a factor of about 1.5 as discussed below (9). In- 
deed, a more restrictive limit of Tr < 10 7 GeV is found 
within the CMSSM by using the full expressions for 
given in [14] which include also the electroweak contribu- 
tions to thermal gravitino production [15, 47, 48]. 

Let us comment on the model dependence of the 
BBN constraints. As described above, the BBN con- 
straints shown in Figs. 1 and 2 are derived with the 



The additional primordial bound-state effects discussed in [43, 
46, 49-53, 94] do not affect the conclusions of this letter; cf. [54]. 
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yield (15). This does introduce a model dependence into 
the Xj^ ax values discussed in the preceding paragraph. 
However, one should stress two points: (i) The yield (15) 
is quite typical for an electrically charged massive ther- 
mal relic [18]. (ii) The -dependent upper limits on Y~- 
are very steep in the relevant region as can be seen, e.g., 
in Fig. 1 of Ref. [47]. Indeed, this steepness is reflected by 
the relatively weak mj dependence of the 6 Li constraints 
that can be seen in Figs. 1 and 2. These two points sup- 
port that the shown BBN constraints and the associated 
ymax vames are q U ite robust. Only for very generous up- 
per limits on the primordial 6 Li abundance and/or excep- 
tionally small Yr- values can the CBBN bound be relaxed 
1 

substantially. The latter could be achieved, for exam- 
ple, with non-standard entropy production after thermal 
freeze out of the h NLSP and before BBN [15, 45, 95] 
which we do not consider since we assume a standard 
cosmological history throughout this letter. 

PROBING T R BY MEASURING THE 
GLUINO-SLEPTON MASS RATIO 

In this section we present a way that could allow us 
to probe with collider experiments the T^ lax value im- 
posed by the 6 Li constraint — or any other constraint that 
can be translated into a similar upper limit on tj- — and 
the dark matter constraint (4). As discussed above, the 
6 Li constraint implies basically an upper limit on the l\ 
NLSP lifetime 

% < r max , (18) 

which is, e.g., r max ~ 5 x 10 3 s for the case of the 
yield (15) as can be seen explicitly in Fig. 2. With a 
given r max and (11) derived from the dark matter con- 
straint (4), one arrives immediately at the following rriQ- 
independent upper limit on the gluino-slepton mass ratio 
at the weak scale: 

/ 2.37 x 10 9 GeV y/ ^pM/i 2 V 

c - V TT X ) V 0.1 J 

f T max \i f m r \* 

x h — = c max , (19) 

VlO 4 sj \100 GcV J ' v ' 

i.e., a reheating temperature of at most T^ ax can viable 
with a given rnj only for a mass ratio c < c max . In other 
words, for a lifetime constraint r max inferred from cos- 
mological considerations (such as CBBN) and/or collider 
experiments, a simultaneous measurement of mj and the 
mass ratio c will provide a model-independent conserva- 
tive limit T™ ax for G LSP scenarios with a l x NLSP. 

In Fig. 3 the solid lines (dashed lines) show c max im- 
posed by the lifetime constraint tj- < 3 x 10 3 s (10 4 s) and 
and the dark matter constraint fl^h 2 < JIdm^ 2 < 0.126 
as a function of mj for values of X^ lax ranging from 




FIG. 3: Upper limits on the mass ratio c = rrig{Mz)/mj i 
imposed by Q. 5 h 2 < n OM h 2 < 0.126 and r Ti < 3x 10 3 s (10 4 s) 
are shown as a function of mj by the solid lines (dashed lines) 
for values of Tg"* ranging from 10 8 GeV up to 3 x 10 9 GeV. 

10 s GeV up to 3 x 10 9 GeV. The lifetime constraint 
tj- < 5 X 10 3 s is slightly weaker (and thereby more con- 
servative) than the 6 Li CBBN constraint obtained with 
the yield (15) for mj^ > 100 GeV as can be seen in 
Fig. 2. Accordingly, the T^ ax values for c = 1 and c = 7 
discussed in the previous section are more restrictive 
than the corresponding T™ ax values in Fig. 3. The even 

more conservative c max curves obtained for rr < 10 4 s 

'i — 

are presented to cover also the case of a more relaxed 
CBBN constraint which could result, for example, from 
6 Li/H| p < 2.7 x 10~ 10 [53], i.e., an upper limit that is 
about an order of magnitude more generous than (16). 
This emphasizes the conservative character of the corre- 
sponding c max curves shown in Fig. 3. 



TESTING THE VIABILITY OF THERMAL 
LEPTOGENESIS AT COLLIDERS 

Thermal leptogenesis provides an attractive explana- 
tion of the baryon asymmetry in the Universe [24] . Since 
successful thermal leptogenesis with hierarchical right- 
handed heavy Major ana neutrinos requires a reheating 
temperature of T R > 10 9 GeV [25-28], one will be able 
to test its viability at the LHC with the method described 
in the previous section. 

From earlier studies of gravitino LSP scenarios [12- 
14, 36], it is known that thermal leptogenesis can be as- 
sociated with testable (mg-dependent) upper limits on 
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the masses of the gluino and the other gauginos. Never- 
theless, the rr constraint from CBBN of 6 Li is not taken 
into account in any of these studies. In fact, as mentioned 
in the Introduction, this rr constraint disfavors thermal 
lcptogenesis in constrained scenarios with universal soft 
SUSY breaking parameters at M G ut [15, 47, 48, 55]. 
Interestingly, in this work, it is this upper limit on tj- 
that allows us to arrive in a model independent way at a 
testable prediction of thermal leptogencsis that does not 
depend on the gravitino mass vHq. 

As can be seen in Fig. 3, already the conservative limits 
imposed by < 10 4 s show that Tg** = 10 9 GeV is as- 
sociated with a gluino-slepton mass ratio of c < 3 in the 
collider-friendly region of < 1 TcV. Thus, this ratio 
favors a gluino mass that will be accessible at the LHC. 
In this way, c < 3 can be understood as a testable pre- 
diction of thermal leptogenesis for phenomenology at the 
LHC. Indeed, the realistic, less conservative constraints 
% < 5 x 10 3 s and T R > 3 x 10 9 GeV point even to the 
exceptional case with a gluino that is only slightly heav- 
ier than the slepton NLSP: c < 1.5. If realized in nature, 
this will guarantee unique signatures at the LHC. 



tj- < 5 x 10 3 s from CBBN of 6 Li implies a truly conser- 
vative limit of Tr < 3 x 10 9 GeV in the collider-friendly 
region of nij^ < 1 TeV and in the limiting case of a 
small gluino-slepton mass splitting at the weak scale: 
c = m-g/mj^ ~ 1. Indeed, with c = 7 — as obtained in 
scenarios with universal soft SUSY breaking parameters 
at the scale of grand unification — our robust conservative 
upper limit is Tr < 10 8 GeV and thus as severe as in sce- 
narios with an unstable gravitino of mass < 5 TeV. 
The rr constraint from CBBN of 6 Li has allowed us 

d 

to derive upper limits on Tr that depend on and c 
only. In particular, we find that the condition for success- 
ful thermal leptogenesis with hierarchical right-handed 
heavy Majorana neutrinos, T R > 10 9 GeV (3 x 10 9 GeV), 
implies an upper limit of c < 3 (1.5) for < 1 TcV. 
This is a prediction of thermal leptogenesis that will be 
testable at the upcoming LHC experiments. 

I am grateful to Josef Pradler for valuable discussions 
and comments on the manuscript. This research was par- 
tially supported by the Cluster of Excellence 'Origin and 
Structure of the Universe.' 



CONCLUSION 

The observation of a (quasi-) stable heavy charged 
slepton li as the lightest Standard Model superpartner 
at future colliders could provide a first hint towards the 
gravitino LSP being the fundamental constituent of dark 
matter. Under the assumption that the gravitino is the 
LSP and a long-lived l\ the NLSP, we have shown that 
measurements of the slepton mass nij^ and its lifetime 
can provide the gravitino mass niQ and upper limits on 
the reheating temperature Tr after inflation, which can 
be tightened once the gluino mass nig is measured. 

The conceivable insights into mg are a direct conse- 
quence of the SUGRA Lagrangian and of R-parity con- 
servation. Thus, they will be valid for any cosmolog- 
ical scenario provided the G is the LSP and a l\ the 
NLSP. The possibility to probe Tr at colliders results 
from the Tr dependence of the thermally produced grav- 
itino density and the fact that this density cannot exceed 
the dark matter density. Thereby, the given upper lim- 
its on Tr rely not only on the SUGRA Lagrangian and 
R-parity conservation but also on the assumptions of a 
standard cosmological history and on (gauge) couplings 
that evolve at high temperatures as expected from their 
standard MSSM renormalization group running. Never- 
theless, under these assumptions, the presented Tr limits 
are very robust since they have been derived in a model 
independent way from a conservative lower limit on the 
(thermally produced) gravitino density. 

Prior to collider measurements of m^, t^, and nig, 
primordial nucleosynthesis can be used to constrain mg 
and Tr. In particular, we have found that the constraint 
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